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Abstract. 

A first principles many-body approach is employed to calculate the band structure 
and optical response of nanometer sized ribbons of SiC. Many-body effects are 
incorporated using the GW approximation, and excitonic effects are included using 
the Bethe-Salpeter equation. Both unpassivated and hydrogen passivated armchair SiC 
nanoribbons are studied. As a consequence of low dimensionality, large quasiparticle 
corrections are seen to the Kohn-Sham energy gaps. In both cases quasiparticle band 
gaps are increased by up to 2 eV, as compared to their Kohn-Sham energy values. 
Inclusion of electron-hole interactions modifies the absorption spectra significantly, 
giving rise to strongly bound excitonic peaks in these systems.The results suggest 
that hydrogen-passivated armchair SiC nanoribbons have the potential to be used 
in optoelectronic devices operating in the UV-Vis region of the spectrum. We also 
compute the formation energies of these nanoribbons as a function of their widths, 
and conclude that hydrogen-saturated ribbons will be much more stable as compared 
to the bare ones. 


PACS numbers: Insert suggested PACS numbers here. 
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Successful synthesis of graphene [1], and discovery of its extraordinary electronic 
properties such as high electron mobility and electrical conductivity [2, 3, 4], has 
tremendously enhanced research interest in the held of low-dimensional systems, in 
general, and graphene-like two-dimensional (2D) systems, in particular. One of the 
drawbacks of graphene is the fact that it is a zero bandgap semiconductor, which limits 
its applications in electronic devices. Therefore, it is of great interest to discover novel 
low-dimensional materials which share the attractive properties of graphene, but, in 
addition, also have a hnite band gap. One such candidate, monolayer of hexagonal 
boron nitride (h-BN) known as ’white graphene’, has been successfully synthesized [5] 
using mechanical exfoliation and chemical vapor deposition techniques, and electrically 
it is a large bandgap insulator. Nanotubes and nanoribbons of boron nitride have been 
synthesized and studied[6, 7, 8, 9]. Carbon and silicon are the two most interesting 
elements in group IV because of their abundance which leads to low costs, and also 
because of the fact that their allotropes possess attractive electronic properties. Like 
carbon, silicon also forms a honeycomb structure called silicene which, not only has been 
studied theoretically [10], but also has been experimentally synthesized in high vacuum 
conditions [11, 12]. Si, unlike C, prefers the sp^ hybridization, therefore, silicene, in 
contrast with graphene, has a buckled bilayer structure. Furthermore, it has not been 
possible to synthesize silicene in free-standing form, thus, limiting its potential for device 
applications. 

Just like carbon, it is well known that the bulk SiC has several allotropes, although 
a layered graphitic one has not been experimentally realized so far. However, theoretical 
calculations have shown that graphitic SiC can be realized by using the technique of 
“extreme hole injection” [13]. Therefore, it is conceivable that it can also be isolated in 
a monolayer form, similar to graphene [14, 15]. Because of the fact that basic unit of 
SiC contains two different atoms, corresponding monolayer will exhibit a hnite band 
gap. Similarly, it may be possible to synthesize SiC-based nanostructures such as 
nanoribbons, nanotubes, and nanodisks etc [14, 15]. From the applications point of 
view, SiC nanostructures will have the advantage that their integration into the present 
Si based technology will be relatively easier. Keeping this in mind, in this work we 
undertake a systematic electronic structure study of SiC nanoribbons (SiCNRs). Next, 
we briehy review the existing studies of SiC, and related structures. 

Consistent with the intuition, various electronic structure studies have predicted 
SiC monolayer to be a semiconductor with a band gap of 2.5 eV[16, 17, 18, 19, 20, 
21, 22, 23, 24, 25]. Although, SiC monolayer has not been realized in the laboratory 
as yet, the ultra thin graphitic SiC, which exhibits light emitting properties, has 
recently been synthesized [26]. SiC micro ribbons [27], crystaline 3C - SiC nanobelts 
[28], bicrystalline SiC nanobelts [29], and crystalline SiC nanoribbons [30] have been 
successfully synthesized. Nanowires of SiC have been synthesized by reacting SiO with 
carbon nanotubes at different temperatures [31]. First principles calculations have shown 
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that hydrogenation of graphitic SiC can lead to formation of stable SiC nanotubes [32], 

As far as theory is concerned, SiC nanoribbons (SiCNRs) have been 
studied extensively using the density-functional theory (DFT) based hrst-principles 
methodology [17, 33, 34, 35, 36, 18, 37, 38, 39, 40, 41], Sun et a/. [17] studied hydrogen 
passivated SiCNRs using Perdew-Wang functional (PW91) under the generalized- 
gradient approximation (GGA) exchange correlation scheme. They found that hydrogen 
passivated armchair SiGNRs (ASiGNRs) are semiconducting with a direct bandgap, 
which oscillates and exhibits a gradual overall increase with the width. For the 
nanoribbons studied, they reported Kohn-sham band gaps in the range 2.3-2.4 eV [17]. 
Sun et a/. [17] also found that unlike their graphene counterparts, hydrogen-passivated 
zigzag SiGNRs (ZSiGNRs) have non degenerate spins, and were near half metallic, with 
the smaller gap being a few meVs, while the larger one being close to 0.6 eV. The ferri- 
magnetic state was reported to be most stable, with an anti parellel spin orientation 
between the two edges [17]. Bekaroglu et a/., [18] based upon their hrst-principles DFT- 
GGA calculations, reported that bare ASiGNRs have indirect band gaps in the range 
1.3-1.7 eV. The same authors[18] as well as Morbec et a/. [41] studied the inhuence of 
vacancies on the electronic structure of SiGNRs using the DFT-GGA approach, and 
found that single (Si and G) vacancies induce magnetism in non-magnetic ASiGNRs, 
and double (Si and G) vacancies induce magnetism in ZSiGNRs. Tang et al. [40] 
studied bare-edge and half bare-edge ZSiGNRs using the DFT-GGA approach and 
reported that the bare-edge ZSiGNRs exhibit magnetic metallic behavior, while half bare 
edge ZSiGNRs are magnetic semiconductors with small band gap. Oxygen and sulfur 
terminated ZSiGNRs were studied by Lopez-Benzanilla et a/. [37] using DFT based local 
spin-density approximation (LSDA), and observed that ZSiGNRs can be tuned from half 
metallic to semiconducting. The inhuence of substitution of edge Si and G atoms with B 
and N atoms in ZSiGNRs was studied by Zheng et a/ [38] using the DFT-GGA approach, 
and they found that such substituted ZSiGNRs exhibit semicounducting behavior. The 
inhuence of doping on the electronic and magnetic properties of SiGNRs have also been 
studied by Gosta et al. [39], who showed that the N doped ZSiGNRs are metallic. 

Generally speaking nanoribbons of various materials show interesting properties 
that are vastly diherent from those of their bulk counterparts due to quantum 
conhnement. The bandgap in armchair nanoribbons (ANRs) oscillates over three 
families according to their width, and approaches a value close to that of the parent 
2D sheet for large widths [42, 17, 43]. Zigzag nanoribbons(ZNRs) of both graphene and 
SiG are spin polarized, and the individual spin gaps can be controlled by application 
of an electric held directed along the width [44, 34, 43]. These unique properties of 
nanoribbons help to tune the band gap on the basis of their widths, and the edge 
geometries. Therefore, a good knowledge of their electronic and optical properties is 
essential for their utilization in opto-electronic devices. 

It is a well known fact that the standard DFT based approaches normally 
underestimate the band gaps of bulk semiconductors because they do not include 
the quantum many-body ehects[45]. However, neglecting many-body ehects in low- 
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dimensional systems such as nanoribbons are likely to cause much more severe errors 
as compared to bulk systems, because electron-correlation effects are enhanced in lower 
dimensions [46, 47, 48, 49, 50]. All the aforementioned electronic structure studies 
of SiCNRs have been based upon DFT with various exchange-correlation functionals 
without accounting for the many-body effects. Furthermore, it will also be of tremendous 
interest to explore the optical absorption in SiCNRs in order to explore the possibilities 
of their applications in opto-electronic devices. With these goals in mind, in this work 
we present a systematic study of ASiCNRs in which we go beyond the mean field DFT 
approach to compute the influence of many-body effects on their band structure within 
the GW approximation. We find significant GW corrections of up to 2 eV to the Kohn- 
Sham band gaps. Furthermore, we also calculate the optical absorption spectra of 
these systems within a Bethe-Salpeter equation (BSE) based approach, which includes 
electron-hole interaction, and thus excitonic effects. The absorption spectra of ASiCNRs 
are found to be dominated by strongly bound exciton peaks with binding energies 
up to 2eV. We also performed calculations of formation energies of these nanoribbons 
for various widths, and conclude that with the increasing widths nanoribbons become 
more stable, and, for a given width, hydrogen-satured ribbons are much more stable as 
compared to the bare ones. 

2. Computational Details 

In earlier works from our group we studied the opto-electronic properties of graphene 
nanoribbons using a semi-empirical Hamiltonian based methodology [51, 52]. In the 
present calculations, we have adopted a plane-wave based first-principles methodology 
as implemented in the ABINIT package. [53, 54] First, the calculations were performed 
for the ground state of the nanoribbon concerned, using DFT-GGA methodology, em¬ 
ploying Perdew, Burke, Enzerhof (PBE) pseudopotentials. [55] During the geometry op¬ 
timization, a kinetic energy cutoff of 26 Hartree (1 Hatree = 27.21 eV) was used, while 
for Brillouin zone integration, a k-point grid with 1x1x11 points was employed. For 
total energy/cell, convergence threshold of 10“^*^ Hartree (Ha) was used, and during the 
structural relaxations, atomic force cutoff of 5 x 10“® Ha/Bohr was employed. Once 
the geometry optimization was complete, charge density profiles, and denisty of states, 
were computed using the software package VASP,[56] employing the same exchange- 
correlation potential, and the kinetic energy cutoff as in ABINIT, but a superior k- 
point mesh of dimensions 1 x 1 x 45. The GGA wave functions are a good starting 
point for the many body perturbation theoretic computations, aimed at including the 
electron-correlation effects. The screening effects are included within a plasmon pole 
approximation, [57] and quasiparticle corrections to the GGA eigenvalues are then eval¬ 
uated using the GOWO approximation. [58] For the exchange and correlation parts of 
self energy calculations, cut off energies of 18 Ha, and 6 Ha, respectively, were used. A 
1 X 1 X 30 k-point grid ensured convergence of the GW band gap to within 0.1 eV for the 
systems considered here. A total of 250 bands was employed for the GW calculations. 
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out of which the number of virtual (unoccupied) bands was > 150. Thus, the choice 
of computational parameters employed in the GW calculations suggests that these cal¬ 
culations are well converged. Lastly, the BSE is solved to obtain the optical response 
with the inclusion of the electron-hole interaction, retaining only the resonant part of 
the Bethe-Salpeter Hamiltonian, because the inclusion of the coupling had a negligible 
effect. In the GW and the BSE calculations, a Goulomb truncation scheme was also 
used to avoid the consideration of long range interactions between the periodic images. 
Nanoribbons were taken to be periodic along the z direction, and to represent an iso¬ 
lated system, vacuum separations of more than 11 A were taken between the adjacent 
unit cells along both x and y directions. 

To denote the widths of the nanoribbons, the usual index notation is used. 
ASiGNRs are prefixed by the number of dimer (bond) lines, W, along the width of 
the ribbon (W-H-ASiGNR is used for H-passivated, and W-ASiGNR used for bare 
armchair nanoribbons). In this work, we have studied ribbons with 4 < W < 12, 
which for H-ASiGNRs corresponds to widths between 0.7 nm and 2.3 nm, and for bare 
ASiGNRs to widths between 0.46 nm and 1.7 nm. 


3. Results and Discussion 

In this section we present the results of our calculations on the quasiparticle band 
structure and the optical absorption spectra of ASiGNRs. Before that, however, we use 
the same methodology to compute the band structure of 2D SiG sheet at the GW level, 
in order to validate our approach. 

3.1. SiC Monolayer 

Because the calculations presented here are the hrst reported GW calculations on 
ASiGNRs, it is important to benchmark the underlying approach. Therefore, we hrst 
use the same methodology to compute the quasi-particle band structure of monolayer 
SiG for which GW calculations had been reported earlier by Louie and coworkers [21]. 
The honeycomb structure adopted for the inhnite SiG sheet is shown in Fig. 1(a), while 
its calculated GW band structure is presented in Fig. 1(b). We obtain a direct band 
gap of 3.4 eV located at the K point of Brillouin zone, at the GW level, which is 0.9 eV 
greater than our DFT-GGA band gap of 2.5 eV. Our GW result is found to be in very 
good quantitative agreement with the GW value of 3.6 eV, also located at the K point, 
reported by Louie group [21]. Furthermore, we also calculated excitonic binding energy 
of 2D SiG sheet using the GW-I-BSE approach. Fig 1 (c) shows the optical absorption 
spectra of sheet, with (red curve), and without (blue curve), electron-hole interactions. 
The difference between the first prominent peaks of these two curves gives the binding 
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(a) (b) (c) 



Figure 1. (Color online) Geometric structures of (a) 2D SiC sheet (d) Hydrogen 
passivated, (e) bare armchair SiC nanoribbons, . Both types of SiCNRs correspond to 
Na = 6, and unit cells employed for the calculations are shown in the dashed boxes. 
Panels (b) and (c) show the calculated GW band structure and optical absorption 
spectra of SiC sheet. Panels (f) and (g) display 2D charge density (in the units of 
e/Bohr^) plots of 6-H-ASiCNR, and 6-ASiCNR, respectively. 


energy of 1.5 eV, which is 0.3 eV more than the value 1.2 eV reported by Louie and 
coworkers [21], This agreement between the GW and BSE calculations performed using 
two different packages, assures us of the quality of our calculations. 

3.2. Quasiparticle energies of hydrogenated armchair SiC nanoribbons 

First we present and discuss the results of hydrogen-passivated ASiCNRs (H-ASiCNRs), 
for which three different families of width, that is, Na = 3p, Na = 3p-|-l, and Na = 3p-|-2 
{p an integer), have been studied, with Na ranging from 4 to 12. The geometrical 
structure of H-ASiCNR for A'^q = 6 is shown in Fig. 1(d), for which our optimized 
average Si-C bond length is 1.78 A, while Si-H and C-H bond lengths are 1.49 A and 
1.09 A, respectively. Table 1 presents the band gaps at the GGA and GW levels. 
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Figure 2. (Color online) For 6-H-ASiCNR: (a) DFT-GGA band structure (pink arrows 
show the charge densities corresponding to VBM and CBM) (b) GW band structures 
(red arrows denote interband transitions), ct/tt character of bands is also pointed out. 
(c) Atom projected density of states (black dashed line denotes the Fermi level) 


and the exciton binding energies for Na= 6, 7, 8. The GW self-energy corrections 
increase their band gaps by 1.6-1.8 eV, when compared to their GGA valnes. Thus, 
the magnitude of the increase in band gaps for H-ASiGNRs due to electron correlation 
effects is signihcantly greater than the 1 eV increase observed in the monolayer discussed 
in the previous section, and for the increase reported for the bulk SiG [21]. The larger 
correction observed for nanoribbons is a consequence of enhanced Goulomb interaction 
effects, and reduced screening, in lower dimensions. 


Next, we specihcally discuss the case of 6-H-ASiGNR (notation Afa-H-ASiGNR), 
because the other two studied H-ASiGNRs exhibit similar behavior. In Figs. 2 (a) 
and 2 (b), we present our calculated GGA and GW quasiparticle band structures of 
6-H-ASiGNR, from where it is clear that the hydrogen passivated ASiGNRs are direct 
band gap semiconductors as their valence band maximum (VBM) and conduction band 
minimum (GBM) occur at the high symmetrical point T in brillioun zone. Partial charge 




















Figure 3. (Color online) (a) Quasiparticle self-energy corrections to the GGA Kohn- 
Sham energies. Corrections to the nearly free electron states are enclosed by ellipses 
(red-dashed), and (b) variation of the GGA band gap and quasiparticle gaps with 
ribbon width for H-ASiCNRs. 

denisty plots for VBM and CBM of 6-H-ASiCNRs are shown in left side panel of Fig 
2 (a). These plots reveal that main contribntion to the VBM charge density is from 
electrons distribnted aronnd the carbon atoms, while electrons from from both Si and 
C atoms contribnte to the CBM charge density. Fig. 2 (c) shows the atom projected 
density of states (PDOS) of 6-H-ASiCNR. PDOS has been plotted both for internal 
atoms, as well as the ones located on the edges. It clearly show that edge and inter¬ 
nal C atoms contribnte mainly to the states states near to the VBM, while internal Si 
atoms contribnte dominantly to form the states close to the CBM. The contribntion of 
hydrogen atoms appears in the condnction band states far away from the CBM. The 
self energy corrections to varions bands are plotted in Fig. 3 (a), and are dependent 
on their natnre. The tt band states associated with the pz orbitals are perpendicnlar 
to the ribbon plane, extending into the vacnnm region. As a resnlt, the many body in¬ 
teractions for the 71 states are less screened, and lead to large qnasiparticle corrections. 
On the other hand, the a band states lie in the plane of the ribbon, and experience 
greater screening. Thns, the qnasiparticle corrections for the a states are smaller than 
the TT states. In H-ASiCNRs, the bands aronnd the Fermi energy comprise only of the vr 
states, which exhibit corrections which are mostly linear as a fnnction of the energy. A 
few points (encircled, cf. Fig. 3a) which deviate from the main linear region correspond 
to the nearly free electron states which are loosely bonnd to the ribbon, and have weak 
many-body interactions resnlting in smaller qnasiparticle corrections [47] . Compared to 
their GGA shapes, the bands not only shift, bnt also stretch by an average 30%-40%, 
and even cross dne to the non-nniform qnasiparticle corrections resnlting from the GW 
calcnlations. The family behavior of the band gaps is preserved at the GW level, as 
the correction is ronghly proportional to the GGA band gap. Fig. 3(b) shows the band 
gap as a fnnction of the ribbon width in GGA and GW approximation for H-ASiGNRs. 
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From Table 1 it is obvious that the family behaviour of DFT and GW band gaps follow 
the heirachy 3p > 3p + 1 > 3p + 2, in good agreement with earlier studies of DFT 
band gaps of H-ASiCNRs,[17] and quite unlike H terminated graphene [46] and boron- 
nitride [49] nanoribbons which follow the hierarchy 3p -f-1 > 3p > 3p -|- 2 


Na 

Egga 

Egw 

AE 

Efe 

6 (3p) 

2.41 

4.18 

1.77 

1.07 

7 (3p + 1) 

2.38 

4.17 

1.79 

1.76 

8 (3p + 2) 

2.31 

3.91 

1.60 

1.93 


Table 1. The values of GGA and quasiparticle band gaps and their difference for 
H-ASiCNRs are tabulated. The last column gives the highest exciton binding energy. 
All energy values are in eV. 


3.3. Optical absorption spectra of H-ASiCNRs 

Based on GW quasiparticle band structures, we calculated optical absorption spectra 
of H-ASiGNRs, with (GW-BSE), and without (GW-RPA) electron-hole interaction, 
and results are presented in Fig. 4. BSE calculations were performed by including 
five valence and five conduction bands. The polarization direction of the radiation is 
assumed to be along the length of the ribbons, as significant quenching occurs in ID 
systems for perpendicular directions. The blue colored curves denote the absorption 
due to inter-band transitions at the GW-RPA level. The red colored curve denotes 
the absorption obtained by solving the Bethe Salpeter equation, including both the 
electron-electron, and electron-hole interactions. 


As compared to the independent particle absorption spectra, BSE based spectra 
exhibit changes both in the positions of the peaks, as well as the line shapes, with 
excitonic states having binding energies in the range 1.07-1.93 eV. We adopt the 
following notation to label interband transitions: Enm denotes the interband transitions 
from the nth valence band to the mth conduction band, with the band indices counted 
from the Fermi level. The underlying band-to-band transitions are also marked with 
arrows in the band structure plot in Fig. 2 (b). The first peak (En) of ribbon arises due 
transitions between the highest valence band, and the lowest conduction band at the 
zone center, while other prominent peak arises due to transitions between the second 
valence and conduction bands. Weak peaks in between include contributions from cross 
transitions between these bands. For the case of Na = 6, 7, E 22 peak is the most intense 
one, while for Na = 8 case, Eu peak acquires signihcant oscillator strength, and becomes 
the most intense peak. The binding energies of the excitons corresponding to the Eu 
peak are 1.07 eV, 1.76 eV, and 1.93 eV, and for Na = 6, 7, and 8, respectively. In 
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Photon energy (eV) 

Figure 4. (Color online) Optical absorption spectra of the H-ASiCNRs, with (red) 
and without (blue) electron-hole interaction. Lorentzian broadening of 0.03 eV has 
been used. 

comparison, our computed value of exciton binding energy of SiC sheet is 1.5 eV, while 
that of bulk 2H-SiC is reported to be 0.1 eV[21]. Thus, exciton binding energy for 
ASiCNRs is comparable to that of the monolayer, but signihcantly larger than that of 
the 3D bulk SiC. This large difference between the exciton binding energies of ASiCNRs 
and bulk SiC can be attributed to their reduced dimensionality, because in reduced 
dimensions, electron-hole interactions get stronger as quantum conhnement increases 
the overlap between electron and hole wave functions, leading to the formation of small 
radius excitons[49]. When we compare the exciton binding energies of nanoribbons 
of different widths, we hnd a strong width dependence, implying that the strength of 
electron-hole interactions is also strongly width dependent. This suggests that one can 
tune the opto-electronic properties of H-ASiCNRs, which are direct band gap materials, 
by manipulating their width, a property highly useful from the point of view of device 
applications. 

3 . 4 . Quasiparticle energies of bare armchair SiC nanoribbons 

In this section, we present and discuss the quasi-particle band structure of bare 
ASiCNRs. Fig 1 (e) displays the relaxed geometric structure of bare 6-ASiCNR, from 
which it is obvious that it undergoes signihcant edge reconstruction because of the 
dangling bonds, and the Si-C bond length at the edge is reduced from 1.78 A to 1.69 
A. The effect of edge reconstruction is also obvious in total charge density plot of the 
ground states of bare 6-ASiCNRs (Fig. 1 (g)), when compared with that of 6-H-ASiCNR 
(Fig. 1 (f)). Furthermore, we hnd that, at all levels of theory, for the same width, the 
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Figure 5. (Color online) (a) DFT-GGA, and (b) GW band structures of bare 6- 
ASiCNR. Red color arrows denote the interband transitions, pink color arrows show 
the VBM and CBM, black color arrows denote the additional valence and conduction 
bands, ct/tt character of bands is also pointed out. (c) Partial density of states for 
6-ASiCNR, Fermi level shown with black dashed lines. 

At the DFT-GGA level, they are smaller by ~0.7 eV , a result in agreement with the 
results of previous theoretical studies, [18] while at the GW level, the difference reduces 
to ~ 0.3 — 0.4 eV. For bare ASiGNRs also we studied the family behavior by considering 
widths of Na= 6 (3p), 7 (3p-|-l), 8 (3p-|-2);, the results of our calculations are presented in 
the Table 2. Figure 5 displays the calculated GGA and GW band structures for Na = 6. 
From the figure and the table it is obvious that, as in the case of hydrogen-passivated 
ribbons, quasiparticle corrections due to the GW approach broaden the band gaps of 
bare ASiGNRs also by about 2 eV. However, unlike the hydrogen-passivated ribbons, 
the bare ASiGNRs exhibit an indirect band gap because the top of the valence band, and 
the bottom of the conduction band, are located at points F and Z of the Brilliouin zone, 
respectively, whereas in H-ASiGNRs both were located at the F point. As shown in Fig. 
5, as compared to H-ASiGNRs, for the case of bare ASiGNRs we can observe that two 
new, relatively flatter bands, appear near both the GBM and the VBM. As far as the 
GBM region is concerned, it is now entirely composed of these two new bands, while for 
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the valence band case, the new bands are slightly away from the VBM. This restructuring 
of the band structure caused by hydrogen removal leads not only to a reduced band gap, 
but also transforms bare ASiCNRs into indirect band gap materials, as compared to H- 
ASiCNRs. Charge density analysis reveals that these additional valence and conduction 
bands are due to the dangling edge bonds, and, because, the electrons involved do not 
participate in the bond formation, they have relatively less interaction with each other, 
leading to flatter bands. Partial charge desity plots for VBM and CBM of bare ASiCNRs 
are presented in Fig 5 (a). A comparison with the corresponding plots of H-ASiCNRs 
reveals that: (a) VBMs of bare ribbons are similar to the case of hydrogen saturated 
ones in that they derive their charge density contributions from the same set of carbon 
atoms, and (b) CBMs of bare ribbons are entirely different from those of the saturated 
ones in that their charged density is almost entirely due to the edge silicon atoms. From 
the PDOS plots of bare ASiCNRs presented in Fig 5 (c), one can deduce that the main 
contribution to the states near VBM is from interior C atoms, while the contribution 
to the states near CBM is mainly due to the edge Si atoms. 



Figure 6. For bare ASiCNRs: (a) quasi-particle self-energy corrections to the 
GGA Kohn-Sham energies. Corrections to the states corresponding to flat bands in 
conduction (valence) band regions are enclosed inside red-dashed circle (ellipse), (b) 
Variation of the GGA and quasi-particle band gaps as a function of the ribbon width. 


In Fig. 6(a), we show the self energy corrections for ASiCNRs, as a function of 
energies. The corrections to the GGA eigenvalues of bare ribbons are qualitatively 
similar to those for H-ASiGNRs, with the corrections to the vr states being larger, as 
compared to those for the a states, and exhibit a linear variation with respect to the 
energy. The states enclosed in red circle, which deviate from the main region, represent 
corrections to the conduction band states localized on the edge atoms. In Fig. 6(b) 
we present the variation of the band gap as a function of the ribbon width, and it is 
clear that bare ASiGNRs follow the same trends as H-ASiGNRs in that the band gap 
exhibits oscillatory behavior. For the largest width considered {Na = 12), GW band 
gap of H-ASiGNR is 0.64 eV larger than that of the SiG sheet, while for bare ribbon the 
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same difference is 0.46 eV. With the increasing width, band gaps of hydrogen passivated 
nanoribbons are expected to approach the band gap of the SiC monolayer, while in case 
of bare ones it is nnlikely, becanse of the presence of flat bands (Fig.5(a)). 


N 

Egga 

Egw 

AE 

Efe 

6 (3p) 

1.74 

3.72 

1.98 

0.62 

7 {3p + 1) 

1.51 

3.48 

1.97 

1.26 

8 (3p + 2) 

1.62 

3.56 

1.94 

2.48 


Table 2. Values of GGA and GW band gaps, and their difference for bare ASiGNRs 
are tabulated. The last column gives the highest exciton binding energy for the given 
width. All energies are in eV units. 


3.5. Optical absorption spectra of bare ASiCNRs 

In Fig. 7 we present the optical absorption spectrnm of bare ASiCNRs, for W = 6 - 8 . 
For the calculations of the BSE spectra, five valence and five conduction bands were 
included. As for the case of H-ASiCNRs, once the electron-hole interaction effects are 
incorporated, the resulting BSE spectrum of bare ribbons looks completely different as 
compared to the band-to-band absorption spectrum computed at the GW-RPA level. 
Similar to the case of H-ASiCNRs, the excitons in bare ribbons are also strongly bound, 
with the binding energies in bare ribbons, for the same ribbon widths, being larger 
than in the hydrogen passivated ones (c/. Tables 1 and 2). We can also find strongly 
bound excitons in bare ASiCNRs with binding energies of in the range 0.62-2.45 eV. The 
reasons behind these large binding energies are the same as in case of H-ASiCNRs, i.e., 
stronger electron-hole interactions due to quantum confinement, and reduced screening 
in lower dimensions. 


One qualitative difference in the optical properties of bare ribbons as compared to 
the hydrogen-passivated ones is that bare ribbons do not absorb at the band gap because 
they exhibit an indirect band gap, as discussed in the previous section. Therefore, the 
prominent peaks displayed in Fig. 7 do not have an En contribution. In all the ribbons 
considered, the first peak corresponds to an Ei^ transition, while the strongest peak is 
due to E 34 transition. Because of the indirect nature of excitons, bare ASiCNRs cannot 
be used in devices such as LEDs, however, they may be useful in transport oriented 
applications. As far as the location of first absorption peak is concerned, both hydrogen- 
passivated and bare nanoribbons follow the pattern E{3p + 2) > E{3p -|- 1) > E{3p). 
Because of the presence of an indirect band gap in bare ASiCNRs, these materials will 
have poor photoluminscence efficiency, and thus will not be useful in light-generation 
applications. 
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Figure 7. Optical absorption spectrum of bare 6-ASiCNRs, with (GW-BSE), and 
without (GW-RPA), electron-hole interactions. Lorentzian broadening of 0.03 eV has 
been used. 


3.6. Energetic Stability 

Energetic stability of hydrogen terminated ASiCNRs is well known from previous 
theoretical studies. [17] In order to check the stability of bare ASiCNRs, we performed 
calculations of their formation energies, and compared them with those of H-ASiCNRs 
. Following formula was employed for computing the edge formation energy per unit 
length of edge of a given nanoribbon. 


Ef — {Ec — NdE2D — 0-5 NhEh2) /2A, 

where E^. is the total energy/cell of the nanoribbon under consideration, E 2 D is the 
total energy/cell of infinite 2D SiC sheet, Eh 2 is the total energy of a hydrogen molecule, 
Nd is the total number of SiC units in the unit cell, is the total number of H atoms 
in the unit cell {Nh = 0, for bare nanoribbons), while L is the length of the edge (in 
A unit). It is easy to see that a lower value of Ef will imply a higher stability of the 
considered nanoribbon. 

Formation energy as a function of nanoribbon width for both bare and H-saturated 
nanoribbons studied here is presented in Fig 8. The fact that both types of nanoribbons 
exhibit little change in the value of edge formation energy with the increasing width, 
implies a weak width dependence. H-ASiCNRs studied here have a much lower average 
formation energy of 0.07 eV/A as compared to 0.76 eV/A of the bare ones, clearly 
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Figure 8. Formation energies of bare ASiCNRs and H-ASiCNRs as a function of 
width 

demonstrating that H-ASiCNRs will be much more stable when compared to the bare 
ones. Therefore, based upon the formation energy data, we conclude that the synthesis 
of H-ASiCNRs in laboratory will be much more favourable as compared to the bare 
ones. 

4. Conclusions 

To summarize, we performed large-scale hrst-principles calculations on armchair-type 
SiC nanoribbons of widths in the range 4 < Af^ < 12, in order to study their electronic 
structure and optical properties. Both hydrogen passivated and bare edge conhgurations 
were considered, sophisticated many-body methods such as the GW approximation 
for the band structure, and BSE for the optical absorption, were adopted, and the 
convergence of the calculations was carefully checked. Therefore, the results presented 
here are fairly accurate, and constitute, to the best of our knowledge, the hrst set of GW 
and BSE calculations for SiGNRs. From our calculations we conclude that hydrogen- 
passivated ribbons are direct band gap semiconductors, and bare edge-ribbons undergo 
signihcant edge reconstruction, and become indirect band gap semiconductors. As 
expected, self-energy corrections widened the band gaps, and excitonic effects modihed 
optical absorption spectra dramatically for both the H-saturated and the bare ribbons. 
Our results predict that quasiparticle band gaps of hydrogen passivated and bare 
ASiGNRs are ~ 2 eV larger than their DFT-GGA band gaps, due to the inclusion of 
many-body effects within the GW approximation. These large quasiparticle corrections 
to the band gaps suggest enhanced Goulomb correlation effects in reduced dimensions. 
Our BSE based calculations of the excitonic effects predict large excitonic binding 
energies in the range of 0.62-2.45 eV in optical absorption spectra of both types of 
ASiGNRs, which is yet another consequence of reduced dimensions. Key results of our 
calculations can be tested in future experiments once these nanoribbons are synthesized 
in laboratories. Theoretically speaking, it will also be of interest to explore other ribbon 
geometries of SiG, and also ribbons with unequal proportions of Si and G. 
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